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A general scheme to calculate dynamical susceptibilities of strongly correlated electron systems within the 
dynamical mean field theory is developed. The approach is based on an expansion over electron hopping around 
the atomic limit (within the diagrammatic technique for site operators: projection and Hubbard ones) in infinite 
dimensions. As an example, the Falicov-Kimball and simplified pseudospin-electron models are considered and 
analytical expressions for the dynamical susceptibilities are obtained. 



In the last decade the main achievements in 
the theory of the strongly correlated electron sys- 
tems are connected with the development of the 
Dynamical Mean Field Theory (DMFT) which is 
exact in the d — oo limit It was shown by 
Metzner and Vollhardt ^,|| that in the d = oo 
limit self-energies are single-site quantities (do 
not depend on wave vector) which leads to a sig- 
nificant simplification. The same is true for the 
four-vertices in the Bethe-Salpeter equation for 
susceptibilities (see G|). 

The aim of this article is to develop a gen- 
eral scheme to calculate dynamical susceptibili- 
ties within a diagrammatic technique for site op- 
erators (projection or Hubbard ones) for strongly 
correlated electron systems described by the gen- 
eral statistical operator 
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where Hq = Hi is a sum of the single-site 
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contributions. Our approach is based on an ex- 
pansion over electron hopping around the atomic 
limit H (see also Q]) instead of an expansion in 
the local interaction 

In the d — oo limit, the lattice problem with 
tfj(r — t') = 5(t — r'ftij/y/d, is mapped onto 
an effective atomic problem with a dynamical 
mean field t?-(r — r') = <% J ct (t — r'). Single-elec- 
tron Green's functions for the lattice G a (us Ul k) = 



^E, (7 1 (lu^) — £fc] and the effective atomic prob- 
lem G { a\uj v ) = [S^ 1 ^) - J a {ijj v j\ 1 are deter- 
mined by the same single-site irreducible parts 
S cr (w I/ ) H and a closed set of equations for 
S tT (w 1/ ) and J a {ijjy) can be written Here, irre- 
ducible parts are contributions to the single-site 
Green's function which cannot be divided into 
parts by cutting one hopping line. 

In a similar way, the expansion for correlation 
functions built on operators A and B in terms of 
the hopping is 
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where thin wavy lines denote electron hopping 
integrals tfJr — r') and arrows denote single- 
particle Green's functions. In eq. (p|), I 



and 



L are single-site quantities, which are the same 
for the lattice and the effective atomic problems 
and are generalized many-particle Green's func- 
tions j| which will be calculated within DMFT. 

To do this, we calculate a two-particle Green's 
function for the effective atomic problem 

(Ta{a 2 ala 4 ) {a) = Q\ 

which, on the other hand, can be written in the 
following way: 
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Combining this equation with the equation for 
the four-vertex of the atomic problem allows us 
to calculate the single-site four-vertex Sim- 
ilarly, we can calculate other single-site vertices 
I I and ^ and, finally, determine the correla- 
tion functions (||) for the lattice. 

As an example, we consider a binary alloy 
model 

Hi = \P+n % - 9 -Pr ni - ±{P+ P-) /in,, (3) 



where Pr- = i ± Sf for the U = pseudospin- 
electron (PE) model, = Ci, P~ = 1 — ej for 
a binary alloy, and ai a — > di, m 



f}fi, Pr = 1- ftft for the Falicov-Kimball (FK) 
model IB . The single-particle Green's function for 
the effective atomic problem is a coherent sum of 
the Green's functions for subspaces S, 
is equal to 
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which allows us to find solutions for S CT (cov) and 

For the model (|) all many-particle Green's 
functions for the effective atomic problem can be 
obtained in an analytical form which allows us to 
calculate dynamical susceptibilities: pseudospin 
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The expression for 0(T, q) coincides with the one 
obtained by Freericks || and it is known that the 
ground state of the model (||) , for a fixed average 
value of the pseudospin, is not uniform and shows 
either commensurate order, incommensurate or- 
der, or phase separation Q. On the other hand, 
in the case of a fixed value of the field h the pos- 
sibility of a uniform first-order phase transition 
(bistability) appears 
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